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With recent developments in simulating nonadiabatic systems to high accuracy, it has become possible to 
determine how much energy is attributed to nuclear quantum effects beyond zero-point energy. In this work 
we calculate the non-relativistic ground-state energies of atomic and molecular systems without the Born- 
Oppenheimer approximation. For this purpose we utilize the fixed-node diffusion Monte Carlo method, in 
which the nodes depend on both the electronic and ionic positions. We report ground-state energies for 
all systems studied, ionization energies for the first-row atoms and atomization energies for the first-row 
hydrides. We find the ionization energies of the atoms to be nearly independent of the Born-Oppenheimer 
approximation, within the accuracy of our results. The atomization energies of molecular systems, however, 
show small effects of the nonadiabatic coupling between electrons and nuclei. 


I. INTRODUCTION 


There have been several recent discoveries 1-5 suggest¬ 
ing that quantum wave functions, which include both 
electronic and ionic degrees of freedom, have many in¬ 
teresting properties that have yet to be explored. This 
includes the development of equations that exactly fac¬ 
torize a wave function into electronic and ionic com¬ 
ponents , 2,6 the disappearance of conical intersections in 
wave functions of model systems , 3 and the use of quan¬ 
tum entanglement to study electronic and ionic density 
matrices . 4 Extending such studies to realistic systems is 
of broad interest and will considerably expand our un¬ 
derstanding of electron-ion systems. However, treatment 
of ab initio electron-ion systems is challenging, and ap¬ 
plications have thus been limited. The most accurate 
simulations of electron-ion wave functions are generally 
done with very specialized wave functions, which are lim¬ 
ited to rather small systems . 7 Methods are also being de¬ 
veloped to treat larger systems with different regimes of 
validity . 8-21 

As a framework to address these problems in general 
realistic systems, we recently demonstrated that quan¬ 
tum Monte Carlo (QMC) can be combined with quan¬ 
tum chemistry techniques to generate electron-ion wave 
functions . 1 We treated realistic molecular systems and 
demonstrated that our method can be scaled to larger 
systems than previously considered while maintaining a 
highly accurate wave function. In the following we ex¬ 
tend our previous work by considering the simulation of 
a larger set of atoms and molecules. We calculate ion¬ 
ization energies and atomization energies that can be di¬ 
rectly compared with previous results for benchmarking 
purposes. 


II. METHOD 

A. Fixed-Node Diffusion Monte Carlo (FN-DMC) 

Diffusion Monte Carlo 22-27 is a projector method that 
evolves a trial wave function in imaginary time and 
projects out the ground-state wave function. For prac¬ 
tical simulations of fermions, the fixed-node approxima¬ 
tion is introduced, which depends only on the set of elec¬ 
tronic positions where a trial wave function is equal to 
zero. This approximation is different than approxima¬ 
tions typically used in quantum chemistry calculations, 
and in this work we demonstrate that we can generate 
high-quality nodal surfaces for a range of systems that 
include full electron-ion wave functions. 

If the trial wave function has the same nodal surface as 
the exact ground-state wave function, FN-DMC will ob¬ 
tain the exact ground-state energy. Approximate nodal 
surfaces can be generated through optimization of the 
full wave function. Such approximate nodal surfaces have 
been tested and validated on a wide range of systems, and 
consistently provide an excellent approximation of the ex¬ 
act ground-state energy, comparable to the state of the 
art in ab initio simulations . 28 In addition, the energies 
generated with FN-DMC are variational with respect to 
the ground-state energy. 

In all but a handful of previous QMC simulations , 29-35 
calculations are performed with nuclei ” clamped” to their 
equilibrium positions. However, such an assumption is 
not fundamentally required by FN-DMC. 


B. Electronic Wave Function and Optimization 

There are several different approaches for generating 
electronic wave functions for a FN-DMC calculation . 36-39 
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Recent advances 38,40,41 have made it possible to simul¬ 
taneously optimize thousands of wave function param¬ 
eters using variational Monte Carlo with clamped nu¬ 
clei. We use an initial guess for the wave function that is 
generated from complete active space self-consistent-field 
(CASSCF) 42 ’ 43 calculations using the quantum chem¬ 
istry package GAMESS-US. 44 The optimized orbitals are 
then used in a configuration interaction singles and dou¬ 
bles (CISD) calculation to generate a series of configura¬ 
tion state functions (CSFs). 45 For the small systems Li + , 
Be 4 ", LiH and BeH, a CASSCF calculation with a large 
active space is used in place of CISD. The multi-CSF 
expansion of the wave function can be expressed in the 
following form, 

iVcSF 

’Fcisd (r;R 0 )= ^ affair; R a ), (1) 

i= 1 

where r refers to the spatial coordinates of all the elec¬ 
trons and R 0 refers to the equilibrium positions of all the 
ions. <j>i{f) and a = {a±,a 2 ,... } are the CSFs and Cl 
coefficients generated from CISD. The cc-pV5Z basis 46 
is used for the atomic systems and the Roos Augmented 
Triple Zeta ANO basis 47 is used for the molecular systems 
except for the smallest system LiH, where the cc-pV5Z 
basis is used. 

After the multi-CSF expansion is generated, we im¬ 
pose the electron-nucleus cusp condition on each molec¬ 
ular orbital 48 and add a Jastrow factor to the wave func¬ 
tion to include electron correlation. 49 Our Jastrow factor 
contains electron-electron, electron-nucleus and electron- 
electron-nucleus terms. The full electronic wave function 
used in FN-DMC is, 

Mr; R) = e^^ciSD (r; R). (2) 

We optimize the CSF and Jastrow coefficients, a and /3, 
respectively, simultaneously with QMCPACK. 50,51 Op¬ 
timization is performed with the ions clamped to their 
equilibrium positions R 0 . The equilibrium geometries 
for BeH and BH are chosen to be the ECG-optimized 
distances for comparison with the ECG (explicitly corre¬ 
lated Gaussian) method, and the geometries for the rest 
of the hydrides are taken from experimental data. 52 We 
use 3.015 a.u. as the equilibrium inter-nuclei distance 
for LiH, as this geometry is found to provide a lower 
clamped-nuclei ground-state energy than the ECG opti¬ 
mized distance of 3.061 a.u.. We include all CSFs with 
coefficients larger than a specific cutoff e to lend reason¬ 
able flexibility to the wave function during optimization. 
We include as many CSFs as possible to maximize the 
flexibility of the wave function. However, the inclusion 
of too many CSFs with small expansion coefficients can 
introduce noise as they require a large number of sam¬ 
ples in the optimization step to be optimized. We have 
chosen e to restrict the number of CSFs in the wave func¬ 
tion to be ^1000 in all systems studied. Optimization 
is performed with the linear method 41 with roughly 10 6 
statistically independent samples. 
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FIG. 1. Dragged-node approach for simulation of atomic and 
molecular systems in QMC. (a) For atomic systems we can 
consider the entire wave function shifting with the ion. This 
process can be visualized by following a contour of the wave 
function. The thick dashed circle represents a contour of the 
electronic wave function when the proton is at its reference 
position R 0 , and the thin dashed circle represents the same 
contour when the proton has moved to a new position R. To 
evaluate the ion-dependent electronic wave function ip e (r, R), 
we simply map the electron to its proper place in the reference 
wave function i/> e (r; R 0 ). That is, ij) e {r,R) = tp e (f+ s,R 0 ) = 
ip e (r + s; R 0 ) where s is the shift required to put the proton 
back to its reference position, (b) For H^", we pick one of 
the protons as an “anchor” and approximate the new wave 
function by dragging the reference wave function with the 
“anchor” proton. We also rotate the wave function to align 
its axis of symmetry with the orientation of the two protons. 


C. Electron-Ion Wave Function 

Once a satisfactory electronic wave function has been 
obtained, we construct the electron-ion wave function us- 
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ing the ansatz, 

Vei(r, R) = il}i{R)^} e {r, R), (3) 

where R denotes the spatial coordinates of all ions and 
ip e (r, R) is an ion-dependent electronic wave function 
adapted from the clamped-nuclei wave function R 0 ) 
through basis set dependence. Due to the localization of 
Gaussian basis sets around nuclei, as used in quantum 
chemistry calculations, the nodes of ^j e change based 
on the ionic positions, which we have previously called 
the dragged-node approximation. 1 Although there are 
approaches for going beyond the dragged-node approx¬ 
imation, it was demonstrated to be highly accurate over 
a range of molecules in previous work. 1 For the systems 
considered here, we can impose various symmetries of the 
Hamiltonian onto the wave function that arise from the 
relative motion of the ions. In Fig. 1 we demonstrate this 
approach for the simple cases of a hydrogen atom and an 
H^~ molecular ion. This approach can be generalized for 
use in larger systems or even applied to parts of a bigger 
system, e.g., treating light ions as quantum particles and 
heavy ions as ’’clamped”. 

The term ifii consists of simple products of Gaussian 
wave functions over each pair of nuclei, 

V’/OR) (X H e -ay(|fl i -^|-6y) 2 ; (4) 

i,j>i 

where a^- is a coefficient that is optimized and bij are 
taken to be the equilibrium distances between the nuclei. 
Since i/ji is nodeless, the choice of the variational param¬ 
eters a.ij and bij does not affect the converged FN-DMC 
energy. FN-DMC is then performed with the fully opti¬ 
mized electron-ion wave function. We perform timestep 
extrapolation for all of the tested systems. At least four 
timesteps from 0.005 Ha -1 to 0.0005 Ha -1 are used for 
all systems studied in the clamped-nuclei FN-DMC cal¬ 
culation, and at least three timesteps from 0.005 Ha^ 1 
to 0.0001 Ha -1 are used in the nonadiabatic FN-DMC 
calculation. 

Using definitions from Ref. 53 , the adiabatic approxi¬ 
mation will refer to the complete neglect of the nona¬ 
diabatic coupling matrix when the Schrodinger equation 
is expressed in the basis of eigenstates of the electronic 
Hamiltonian. In this context, the nonadiabatic contri¬ 
bution to an eigenvalue of the electronic Hamiltonian 
can be partitioned into two parts: the diagonal Born- 
Oppenheimer correction (DBOC), which only involves 
the single electronic state of interest, and the remain¬ 
ing corrections arising from terms that involve excited 
eigenstates of the electronic Hamiltonian. The DBOC 
discussed in this work is the expectation value of the 
nuclear kinetic energy operator for the ground adiabatic 
electronic state. We define the clamped-nuclei ground- 
state energy E c as the lowest eigenvalue of the electronic 
Hamiltonian and the nonadiabatic ground-state energy 
E n as the lowest eigenvalue of the full molecular Hamil¬ 
tonian that includes the nuclear kinetic energy. The zero- 
point energy (ZPE) for a diatomic molecule is the energy 


of the ground vibrational state of the one-dimensional vi¬ 
brational mode. Note that the ZPE of the nuclei is part 
of the difference E n — E c . The ZPE is not considered to 
be nonadiabatic, but its contribution is included in the 
full molecular Hamiltonian. 


III. RESULTS AND DISCUSSION 
A. Atoms and Ions 

To assess the quality of our results for atoms and 
ions 67 , we compare to previous results from highly ac¬ 
curate simulations, as presented in Table I. For the 
clamped-ion results, QMC 37 ~ 39, 68,69 and quantum chem¬ 
istry benchmarks are available for comparison. To illus¬ 
trate the high-quality QMC techniques used in this work, 
we compare our clamped-ion atomic results with a recent 
QMC benchmark study. 39 The ground-state FN-DMC 
energies consistently agree across all systems studied (ex¬ 
cept for 0+) within 0.1 mHa. This shows that similar 
nodes can be obtained with different forms of the wave 
function. In particular, our large (~ 1000 CSF) multi¬ 
determinant expansions can be compared with the ap¬ 
proach used by Seth et al ., 39 which relies on moderately- 
sized multi-determinant expansions (~ 100 CSF) with a 
backflow transformation. For certain atoms we can com¬ 
pare to more accurate simulation techniques. For C + 
as well as the neutral and ionized Li, Be and B, highly 
accurate ECG calculations that are all converged well be¬ 
yond 0.1 mHa to the true ground-state energy are avail¬ 
able. The convergence is corroborated by results from 
the Hylleraas method for Li 54 and Be + . 62 In Table I we 
have used the lowest variational results as our references 
for these systems, as the convergence is such that the 
accuracy is higher than other current theoretical or ex¬ 
perimental estimates. 

All of our clamped-ion results agree within 0.2 mHa of 
the ECG references, as shown in Figure 2. The error bars 
for the reference ECG results are absorbed into the DMC 
error bars for clarity, although the ECG error bars are 
orders of magnitude smaller compared to the DMC error 
bars. While ECG results exist for C and N, they are not 
well converged and are not suitable references. 60,70 The 
benchmark results in Ref. 57 are a standard for atomic en¬ 
ergies, and we report them as our references in Table I for 
the larger atoms. However, these benchmark results are 
not consistently accurate to 0.1 mHa. For instance, if we 
use the ECG results for C + with the most accurate ion¬ 
ization reference energy, then we find a reference energy 
for the C atom of -37.84489 Ha, which is 0.1 mHa higher 
than that reported in Ref. 57 . The systems with the most 
error are O and F, for which other QMC studies seem to 
experience similar difficulties. 38,39,71,72 We note that for 
some of these systems it may be possible to absorb the 
sign problem and increase the accuracy further in future 
studies. 73,74 

It is more difficult to find accurate references for the 
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TABLE I. Ground-state energies for atoms and ions and the ionization energies for the atoms: fixed-node DMC 
results of this work (FN-DMC) for atoms and ions with and without the Born-Oppenheimer approximation. The 
rows marked with bold FN-DMC are our nonadiabatic results. The ionization potentials (IPs) are reported in 
the last section of the table. Energies are given in units of Hartree. For the highly accurate Hylleraas and ECG 
results, up to 8 digits are reported in the table. 

Atom 

Li( 2 S) Be( b S) B( 2 P) 

C( 3 P) N( 4 S) 0( 3 P) F( 2 P) 

FN-DMC 

Seth DMC 39 

tr a 
ref 

FN-DMC 

ECG b 

clamped-ion 

-7.478057(5) -14.66731(1) -24.65374(2) 

-7.478067(5) -14.667306(7) -24.65379(3) 
-7.4780603 -14.667356 -24.653866 

nonadiabatic 

-7.47742(2) -14.66643(3) -24.65252(4) 

-7.4774519 -14.666435 -24.652624 

-37.84448(2) -54.58851(6) -75.0658(2) -99.73177(6) 
-37.84446(6) -54.58867(8) -75.0654(1) -99.7318(1) 
-37.8450 -54.5892 -75.0673 -99.7339 

-37.84273(4) -54.58641(5) -75.06313(6) -99.7293(1) 
-37.841621 N/A N/A N/A 

Ion 

Li+( 1 S) Be+( 2 S) B+( 1 S) 

C+( 2 P) N+( 4 S) o+( 3 p) F+( 2 P) 

FN-DMC 

Seth DMC 39 

£ref C 

clamped-ion 

-7.27989(2) -14.324749(7) -24.34883(1) 
-7.279914(3) -14.324761(3) -24.34887(2) 
-7.2799134 -14.324763 -24.348884 

-37.43071(2) -54.05371(5) -74.56597(6) -99.0909(1) 
-37.43073(4) -54.05383(7) -74.56662(7) -99.0911(2) 
-37.430880 -54.0546 -74.5668 -99.0928 


nonadiabatic 


FN-DMC -7.27931(4) -14.32387(2) -24.34758(3) -37.42899(6) -54.05165(4) -74.5634(1) -99.0885(1) 

ECG d _N/A -14.323863 -24.347641 -37.429169 N/A_N/A_N/A 


clamped-ion 


IP (FN-DMC) 

0.19817(2) 

0.342 56(1) 0.304 90(2) 

nonadiabatic 

0.413 77(3) 

0.534 79(8) 

0.4998(2) 

0.6409(1) 

IP (FN-DMC) 

0.19811(4) 

0.342 57(4) 

0.304 94(5) 

0.413 74(7) 

0.534 76(7) 

0.4998(1) 

0.6408(1) 

IP (Ref .) 6 

0.198 130 

0.342 572 

0.304 980 

0.414 014 

0.534 775 

0.500 452 

0.640 946 


a For the atomic references, we use the Hylleraas result for Li, 54 and ECG results for Be JJ and B.° 6 Ref. 57 is used 
for C,N,0 and F where the ground-state energies are taken from Table XI. 
b We use nonadiabatic ECG results as the reference for Li, 58 Be 59 and B 56 , which are converged to the true 
ground-state to well within 0.1 mHa. The result for C, 60 however, may have error on the order of 1 mHa. 
c For the ionic references, we use the ICI result for Li + , 61 Hylleraas result for Be" 1 " 62 and ECG results for B +63 
and C+. 7 ' 64 Ref. 57 is used for N+,0+,F + . 
d ECG references only exist for Be + , 59 B +63 and C + . 64 

6 Spin-orbit coupling and relativistic corrections 65 are removed from experimental data 66 for comparison. 


nonadiabatic results. We provide the first nonadiabatic 
QMC benchmarks for the first-row atoms. There are six 
ECG calculations of nonadiabatic ground-state energies 
that are reportedly converged beyond 0.1 mHa, which 
we use as references. Our reported nonadiabatic ground- 
state energies for Li, Be, Be + , B, B + and C + are in 
agreement with the ECG results to within 0.2 mHa, as 
shown in Figure 2. For these systems, the ECG results 
are converged to essentially the exact ground-state en¬ 
ergies in both the clamped-ion and nonadiabatic cases. 
The difference between our DMC ground-state and ECG 
reference is the fixed-node error present in our wave func¬ 
tions. We would expect the clamped ion results to be 
more accurate than the nonadiabatic results, since the 
nonadiabatic wave functions are inherently more difficult 
to construct. However, for the systems in Figure 2, this 
difference in quality is less than 0.1 mHa. In the case of 
Be, Be + , and B, the nonadiabatic wave function is actu¬ 
ally more accurate than the corresponding clamped-ion 
wave function. 

No reference calculations exist for the heavier atoms 
N,0, and F. However, it is possible to apply finite-mass 


correction 57 ’ 75 (i.e., divide by l + m e /M, where m e is the 
mass of an electron and M is the mass of the nucleus) to 
the best clamped-ion references to estimate the nonadi¬ 
abatic references. The energies for N, O, and F obtained 
in this way are -54.5871, -75.0647 and -99.7310 Ha, re¬ 
spectively. For the ionized states, we obtain -54.0525, 
-74.5643 and -99.0900 Ha. 

The ionization potentials are reported in Table I and 
shown in Figure 3. For determining a set of nonadiabatic 
reference data, we subtract the spin-orbit and relativis¬ 
tic corrections (estimated by Klopper et. al. 65 ) from the 
NIST experimental data. 66 Ref. 65 is considered to have 
the most accurate ionization energies due to its usage of 
state-of-the-art quantum chemistry techniques shown to 
provide close agreement with experiment. For the atoms 
considered in this work, ionization energies have previ¬ 
ously been predicted to be independent of all nonadia¬ 
batic effects beyond the DBOC to within an accuracy of 
0.1 mHa. 65 This prediction is based on calculations that 
are reported to be exact and agree to high accuracy with 
experiment. As shown in Figure 3, the ionization poten¬ 
tials calculated with and without the Born-Oppenheimer 
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FIG. 2. FN-DMC ground state energies for Be + , Be, B + , B, 
C + relative to ECG references 55,56,59,62-64 for either clamped- 
ion or nonadiabatic calculations. These relative energies pro¬ 
vide an estimate for the fixed-node error in the electronic and 
electron-ion wave functions, respectively. 



FIG. 3. Calculated ionization energies relative to reference 
data. The same reference is used for both clamped-ion and 
nonadiabatic results. The calculated energies are all within 1 
mHa of the reference. 


approximation are all within 1 mHa of the reference ener¬ 
gies. Further, the clamped-ion and nonadiabatic predic¬ 
tions for the ionization potentials are statistically indis¬ 
tinguishable for all systems studied, consistent with the 
previous study. 65 

In Table II and Figure 4, we demonstrate the amount 
of nonadiabatic contribution to the ground-state ener¬ 
gies in atoms and ions calculated as the difference be¬ 
tween the nonadiabatic and clamped-ion ground-state 
energies. The amount of nonadiabatic contribution is al¬ 
ways positive for these systems and mostly increases with 
atomic number. Using previous benchmark values for the 
DBOC, we can break down the nonadiabatic contribution 
of our system into a DBOC contribution and everything 
beyond the DBOC. 77-79 The DBOC is relatively insen¬ 



FIG. 4. The nonadiabatic contribution to ground-state ener¬ 
gies of atoms and ions calculated with FN-DMC. The nona¬ 
diabatic contribution is partitioned into the DBOC and the 
remaining correction. A hatched bar indicates the contribu¬ 
tion is negative. The numerical DBOC data is provided in 
Table II. 



FIG. 5. The nonadiabatic contribution to ground-state en¬ 
ergies of ions and their corresponding hydrogen-like atoms 
calculated with FN-DMC and analytically as shown in Eq. 5. 


sitive to the level of theory. Figure 4 indicates that in 
the atomic systems, the DBOC is the dominant contri¬ 
bution to the nonadiabatic energy, with the remaining 
amount being close to zero within error bars. The nona¬ 
diabatic energy is relatively constant between the neutral 
and cationic species. This observation suggests that the 
amount of nonadiabatic contribution is insensitive to the 
addition or removal of a valence electron. Physically, the 
valence electrons are farther from the nucleus than the 
core electrons, and thus are likely to be affected to a lesser 
degree by the delocalization of the nucleus. 

The nonadiabatic contributions in the cations can also 
be compared with those in their corresponding hydrogen- 
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TABLE II. Nonadiabatic corrections for the ground-state en¬ 
ergies of atoms and ions. E n and E c are the FN-DMC calcula¬ 
tions of the nonadiabatic and clamped ground-state energies, 
respectively. The DBOC contribution is provided by Wim 
Klopper.' 6 All energies are reported in units of mHa. 


System E n — E c 

DBOC 

System E n — E c 

DBOC 

Li+ 

0.58(4) 

0.591970 

Li 

0.64(2) 

0.608411 

Be+ 

0.88(2) 

0.899706 

Be 

0.88(3) 

0.920848 

B+ 

1.25(4) 

1.242988 

B 

1.21(5) 

1.241669 

C+ 

1.72(6) 

1.710382 

C 

1.75(5) 

1.710900 

N+ 

2.07(6) 

2.066914 

N 

2.10(8) 

2.069149 

0+ 

2.6(1) 

2.440320 

O 

2.6(2) 

2.441821 

P+ 

2.4(2) 

2.675128 

F 

2.5(1) 

2.678181 


like atoms for a more in-depth analysis. The nonadia¬ 
batic contribution in a hydrogen-like atom can be ob¬ 
tained analytically. The result in Hartree atomic units 
is 

E n — E c = “2~(1 — A 4 ) (5) 

where fi = is the reduced mass of the hydrogen-like 
atom and M and Z are the mass and atomic number 
of the nucleus, respectively. The increase in the nonadi¬ 
abatic contribution with increasing Z for hydrogen-like 
atoms reflects the stronger Coulombic attraction between 
the electron and the nucleus, which enhances the effects 
of the delocalization of the nucleus. An interesting case 
to consider is the transition from Li 2+ to Li. As shown in 
Figure 4 and Figure 5, the addition of a core electron to 
Li 2+ decreases the nonadiabatic contribution, while the 
addition of a valence electron has no further effect within 
our error bars. We also calculate the nonadiabatic contri¬ 
bution in Be 2+ to be 0.78(5) mHa, which is 0.29(5) mHa 
lower than the nonadiabatic contribution in Be 3+ and is 
closer to that in Be + of 0.88(2) mHa. Because the core 
electrons interact more strongly with the nucleus than do 
the valence electrons, the core electrons are affected more 
by the delocalization of the nucleus. Moreover, the addi¬ 
tion of a second core electron decreases the nonadiabatic 
contribution for Li 2+ and Be 3+ . We note that the nona¬ 
diabatic correction to the atomic ground-state energies 
of Eq. (5), which only holds for single electron systems, 
is roughly linear in Z, while the relativistic recoil correc¬ 
tion 80 scales as Z 4 . Therefore, the nonadiabatic effect is 
not seen experimentally, as it is less significant than this 
relativistic effect. 


B. Hydrides 

In Table III, we present our results on a series of molec¬ 
ular systems (hydrides). Finding accurate reference data 
for these systems to 0.1 mHa is not straightforward. We 
will use highly converged ECG data when available. Two 
ECG calculations have been performed in the clamped- 
nuclei limit for LiH 75,84 and we agree within 0.03 mHa 


with the more recent reference. For the rest of the sys¬ 
tems, we combined the best clamped-ion atomic refer¬ 
ences in Table I and thermochemistry 85 estimates of at¬ 
omization energy D e in Table III to produce the reference 
ground-state energies. For BeH and BH, we are within 1 
mHa of the reference values, and our energies are lower 
than the best available quantum chemistry results of - 
15.247846 Ha 86 and -25.287650 Ha 87 for BeH and BH, 
respectively. 
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FIG. 6. The nonadiabatic FN-DMC ground-state energies of 
LiH, BeH and BH relative to ECG references. The error bars 
for the nonadiabatic ECG references are shown as thick dark 
lines, and the error bars for the FN-DMC calculations are 
comparable to the size of the symbols. 

Nonadiabatic ECG calculations only exist for the three 
smallest hydrides. Our results for LiH and BeH agree 
with the ECG references to within 0.2 mHa, as shown 
in Figure 6. The ECG reference for LiH is converged to 
the true ground-state energy beyond 0.1 mHa; thus, it is 
likely that our wave function has a fixed-node error of 0.2 
mHa. For BeH, our result is within 0.1 mHa of the ECG 
reference and agrees within error bars. With BH being 
one of the largest ECG simulations performed, the DMC 
result is actually lower in energy, in this case by 1 mHa. 
The ECG error bar on BH is large, and it is not evident 
how close our result is to the true ground state, although 
extrapolating the ECG result with basis set size suggests 
we are within 1 mHa. 59 For these nonadiabatic systems, 
we have the lowest variational result for BH, and the only 
simulated results of for CH, OH, and HF, to the best of 
our knowledge. 

The atomization energies of the diatomic systems 
are reported in Table III. High-quality thermochemistry 
benchmarks are used for comparison. 85 We take the refer¬ 
ence energies from the last column of Table VI of Ref. 85 
and subtract the corrections in the A Esr (scalar rel¬ 
ativistic) and SO (spin-orbit coupling) columns for the 
comparison with our nonadiabatic energies. For the com¬ 
parison with our clamped-nuclei results, we further sub¬ 
tract the DBOC and ZPE (zero-point energy) correc¬ 
tions. Corrections from spin-orbit coupling and relativis- 


X - « 

"T*. 


i v 
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LiH BeH BH 
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TABLE III. Ground-state energies and atomization energies: fixed-node DMC results of this work for 
all first row hydrides with and without the Born-Oppenheimer approximation. The rows marked with 
bold FN-DMC are our nonadiabatic results. All atomization energies are estimated for OK. D 0 includes 
zero-point energy contribution, while D e does not. Both total energies and dissociation energies are 
given in units of Hartree. 


Molecule 

LiH( 1 E+) BeH( 2 E+) 

BH( 4 E + ) 

CH( 2 n) 

OH( 2 n) 

HF( 4 E+) 


clamped-nuclei 




FN-DMC 

rp a 

ref 

-8.070518(7) -15.24793(2) ■ 

-8.0705473 -15.2483(4) 

-25.28867(3) 

-25.2893(2) 

nonadiabatic 

-38.4780(1) 

-38.4792(2) 

-75.7356(1) 

-75.7382(2) 

-100.4552(1) 

-100.4600(3) 

FN-DMC 

ECG 59,81,82 

-8.06624(3) -15.24194(5) ■ 

-8.0664371(15) -15.24203(10) ■ 

-25.28128(9) 

-25.2803(10) 

-38.4672(3) 

N/A 

-75.7245(5) 

N/A 

-100.4431(4) 

N/A 

D e (FN-DMC) 
D e Feller b 

t 

0.092 46(1) 0.080 62(2) 

0.092 62(5) 0.0809(4) 

clamped-nuclei 
0.134 93(3) 
0.1354(2) 
nonadiabatic 

0.1335(1) 

0.1342(2) 

0.1699(2) 

0.1709(2) 

0.2234(1) 

0.2258(3) 

D 0 (FN-DMC) 

D 0 Feller c 

D 0 Exp. 52 ' 83 

0.08910(4) 0.075 78(6) 

0.089 40(5) 0.0761(4) 

0.088 74(38) 0.074 75(4)' 

0.1290(1) 

0.1299(2) 

3 0.1281(37) 

0.1248(3) 

0.1276(2) 

0.1275(5) 

0.1617(5) 

0.1622(2) 

0.1622(1) 

0.2141(4) 

0.2166(3) 

0.2158(3) 


a For LiH, ECG provides the best reference energy. 84 For the rest of the systems, we combined the 
best clamped-ion atomic references in Table I and thermochemistry estimates of D e in this table to 
produce the reference ground-state energies. 

b Estimates for D e are calculated by subtracting the scalar relativistic, spin-orbit coupling and zero- 
point energy corrections from the reference D 0 in Table VI of Ref. 85 . 
c Here only the scalar relativistic and spin-orbit coupling corrections are subtracted. 
d The atomization energy for BeH in Ref. 52 disagrees with previous high-level theoretical bench¬ 
marks, 59 ’ 80 thus we use Ref. 83 instead. For several of the systems, multiple experimental values 
are available in the literature. We report experimental values that were aggregated in one single 
reference, 52 except for BeH. 83 


tic effects are not used, as they are not included in our 
Hamiltonian. The atomization energies estimated in the 
clamped-nuclei limit agree within 1 mHa of the references 
for all but the largest molecule, HF. Within quantum 
Monte Carlo, it is generally more difficult to obtain an 
accurate nodal surface for a molecule than for an atom. 
As a result, our estimates for the clamped-nuclei atomiza¬ 
tion energies are lower than the references in all cases. A 
similar trend can be observed when comparing our nona¬ 
diabatic results with the references. For each molecule, 
the deviation from the reference is similar in the clamped- 
nuclei and nonadiabatic cases except for CH. 

In Figure 7, we compare both our clamped-nuclei and 
our nonadiabatic results to experimental data. We also 
provide adiabatic estimates by adding the zero-point 
energies calculated with coupled-cluster techniques in 
Ref. 85 to our clamped-nuclei results. To calculate exper¬ 
imental atomization energies starting from the clamped- 
nuclei results, energetic corrections due to zero-point mo¬ 
tion of the nuclei, nonadiabatic effects, spin-orbit cou¬ 
pling and relativistic effects should be included. For these 
highly adiabatic systems, the inclusion of zero-point mo¬ 
tion alone is sufficient to bring our clamped-nuclei results 
to within 2 mHa of the experimental results. Except for 
the case of CH, the nonadiabatic results agree closely 
with their adiabatic counterparts and are closer to the 
experimental values, although for BH the experimental 



FIG. 7. Atomization energies of first row hydrides obtained 
with FN-DMC relative to experimental data. The adia¬ 
batic results are estimated by adding zero-point energies from 
Ref. 85 to the clamped-nuclei results. 


error bar is too large to provide a high-accuracy com¬ 
parison. For CH, the experimental result suggests that 
our electron-ion wave function for this molecule has an 
unusually large fixed-node error. 

To estimate the nonadiabatic contribution to the 

















TABLE IV. Nonadiabatic corrections for the ground-state en¬ 
ergies of diatomic molecules. E n and E c are the FN-DMC 
calculations of the nonadiabatic and clamped ground-state 
energies, respectively. The ZPE and DBOC contributions are 
provided by David Feller. 91 The nonadiabatic correction for 
the dissociation energy estimated with FM-DMC are included 
in the A D 0 column. All energies are reported in units of mHa. 


System E n — E c 

ZPE 

DBOC A D 0 

LiH 

4.28(3) 

3.17 

0.902410 -0.19(4) 

BeH 

5.99(6) 

4.65 

1.251000 -0.19(6) 

BH 

7.39(9) 

5.34 

1.692559 -0.6(1) 

CH 

10.8(3) 

6.44 

2.109487 -2.3(3) 

OH 

11.1(5) 

8.43 

2.670397 0.2(5) 

HF 

12.0(4) 

9.34 

2.799624 0.1(4) 


FIG. 8. The nonadiabatic contribution to the ground-state 
energies in hydrides calculated with FN-DMC. The adiabatic 
reference energies are calculated by adding zero-point energy 
contributions from Ref. 85 to our clamped-nuclei results. The 
nonadiabatic contribution is partitioned into the DBOC and 
the remaining correction. A hatched bar indicates the contri¬ 
bution is negative. 


ground-state energies for these hydrides, we calculate 
the difference between our nonadiabatic and adiabatic 
results, as shown in Figure 8. Similar to the atomic 
case, we break down the nonadiabatic energy of our sys¬ 
tem into a DBOC contribution and everything beyond 
the DBOC. 78 - 79 ’ 88 The ZPE and DBOC contributions to 
this difference arc listed in Table IV. We also calculate 
the nonadiabatic correction to the dissociation energies 
of the hydrides. For BeH, OH, and HF, the nonadia¬ 
batic contribution is almost entirely accounted for by the 
DBOC with the remaining correction being zero within 
error bars. For LiH, BH, and CH, the remaining amount 
of nonadiabatic contribution seems to be nonzero, and 
appears quite significant in CH. However, if the electron- 
ion wave function is significantly lower in quality than 
the electronic wave function for a given system, then 
the amount of nonadiabatic contribution will be overesti¬ 
mated. We also use the zero-point energies from Feller et. 
al. as corrections, which may introduce some additional 
uncertainty. Regardless, our current predictions suggest 
that nonadiabatic effects in BH and CH are larger than 
in the other systems we considered. 

For the LiH molecule, we also calculated the electron 
affinity for comparison to ECG results. We calculated 
the ground-state energy of LiH - to be —8.08222(2) Ha 
for the case of clamped-nuclei. With nonadiabatic effects 
included, our result is —8.07811(3) Ha. Our nonadiabatic 
result is in good agreement with a previous ECG study, 89 
which reported a value of —8.07856887 Ha. We report an 
electron affinity of 0.01187(4) Ha, which can be compared 
to the ECG prediction of 0.012132(2) Ha and agrees with 
the experimental value of 0.0126(4) Ha. 90 


IV. CONCLUSION 

We calculated the ground-state energies of first-row 
atoms and their corresponding ions and hydrides with 
and without the Born-Oppenheimer approximation. In 
addition, we examined the amount of nonadiabatic con¬ 
tribution to the ground-state energies of all systems stud¬ 
ied and determined the amount to be up to a few mHa. 
In the case of CH, the nonadiabatic effects beyond the 
DBOC were unusually large, although we could not rule 
out the possibility that this discrepancy is due to the 
fixed-node error in our simulations. We found the ion¬ 
ization energies of the atoms to be independent of the 
Born-Oppenheimer approximation, consistent with a pre¬ 
vious high-level quantum chemistry study. 65 In contrast, 
the atomization energies of the hydrides showed effects of 
nonadiabaticity, although they were generally much less 
than 1 mHa. This work obtained the first nonadiabatic 
QMC benchmark data for non-relativistic ground-state 
energies and obtained the lowest variational result for 
BH and the only results for CH, OH and HF, to the best 
of our knowledge. 

In comparing to accurate benchmark results obtained 
with other methods, we have demonstrated the validity of 
our wave function ansatz, namely it does produce a high- 
quality electron-ion wave function. This technique also 
has the potential to solve interesting larger-scale prob¬ 
lems due to its ease of implementation, as well as the 
polynomial scaling in computational time with respect 
to the number of electrons. 
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